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Abstract: We present the results of a new Hubbard model for optimally doped 
La2-xSrxCu04^. This model uses parameters derived from BLYP calculations on the cluster 
CuOq. It explicitly includes the Cu dx2_y2 and orbitals, the O orbitals, and the apical 
O Pz orbitals. When correlation is properly included in the Hubbard model, we find that 
there is a crossing of two bands in the vicinity of the Fermi level for the optimally doped 
superconductor. This crossing rigorously occurs along the (0, 0) — (vr/a, vr/a) direction of the 
2-D Brillouin zone. The crossing arises due to the overlap of a broad "-Big" band dominated 
by Cu dx'2-y2 character and a narrower "^ig" band dominated by Cu d^^ character. We 
conclude that optimal doping of La2-xSrxCuO^ and related materials is achieved when the 
Fermi level coincides with this crossing. At this point, formation of Cooper pairs between 
the two bands (i.e. inter-band pairing or IBP) leads to superconductivity. Furthermore, 
using geometric considerations, we extend our conclusions to Y Ba2Cu'^0Q^s and offer a 
simple explanation for the seemingly complex behavior of as a function of doping in this 
material. This behavior can be understood on the basis of multiple band crossings. 
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I. Introduction. 

Eleven years have elapsed since the discovery of high-temperature copper oxide 
superconductors.^ In this work, we focus on the recent suggestion by Tahir-Kheli^ that 
superconductivity arises in these materials due to an inter-band pairing mechanism (IBP). 
That is, Cooper pairing of electrons belonging to two distinct bands may be the cause of 
superconductivity. Such a proposal is appealing as it readily explains why superconductivity 
is only observed at very specific doping levels. La2-xSrxCu04^, for instance, shows an opti- 
mal Tc of 39 K when x = 0.15, but superconductivity quickly vanishes as the doping level is 
changed.'^ The IBP theory only requires that two bands cross. Should such a crossing exist, 
then optimal doping is achieved when the Fermi level coincides with this crossing. Precise 
doping levels are needed to achieve this. 

While such band crossings have not been previously noted by others, we find a cross- 
ing to occur in optimally doped Lai g^SrQ i^CuO^ within 0.15 eV of the Fermi level using a 
simple Hubbard model that includes the effect of electron correlation. The crossing occurs 
between one band which is dominated by Cu dx2_y2 character and a second which is domi- 
nated by Cu dg2 character. Moreover, we find that when the model is empirically adjusted 
to include missing electronic effects, the crossing can be seen to occur at exactly the Fermi 
level. As detailed in the accompanying article by Tahir-Kheli,^ a number of key experimen- 
tal observations that arc otherwise anomalous, such as the temperature dependence of the 
NMR spin relaxation rates and Knight shifts, the Hall effect, the resistivity, and Josephson 
tunneling, are easily explained with the resulting band structure and the IBP model. 

What follows is a detailed account of how this Hubbard model was developed for 
Lai_s5Sro_i5Cu04. The resulting band structure is discussed in length. In addition we 
find that other copper oxide superconducting materials can be understood based largely on 
geometric constraints. As an example, we offer simple arguments as to why YBa2CusOQ^s 
appears to show two characteristic T^'s over an extended doping range. 

II. Calculational Details. 
II. A. Cluster Calculations 

Parameters for the Hubbard model were extracted from restricted open-shell density 
functional calculations on the cluster CuOq. These calculations used the gradient corrected 
BLYP functional^ with the standard 6-31+G* basis set^ on the oxygen atoms and Hay 
and Wadt's^ effective core potential and basis set on the copper. All calculations were 
performed using the Jaguar^ ab initio electronic structure program on a dual processor 200 
MHz Pentium Pro running Linux. 

The CuOq cluster was embedded in a point charge array of 1364 ions. The ions had 
the formal charges of +2.000 for Cu, -1.925 for O, -2.000 for apical O, and +2.925 for La/ Sr. 
The total cluster and point charge array had the symmetry of the tetragonal unit cell 
and was 5 unit cells wide (18.940 A) in the a and b directions and 3 unit cells tall (39.618 
A) in the c direction. Fractional charges were used at the edges. It should be noted that at 
the low temperatures where superconductivity appears, the crystal shows a C2/1 distortion 
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to an orthorhombic unit cell. This should be a relatively small perturbation, so for the 
sake of simplicity the higher symmetry structure was used. However, we note here that the 
distortion will have some important implications in regard to superconductivity. This will 
be discussed below. The tetragonal crystal structure was taken from Hazen^ and is given in 
Table I. 

To obtain the Hubbard parameters, the density of a single state was optimized and 
the resulting orbitals were localized. By then making specific combinations of the localized 
orbitals and not allowing them to relax in subsequent DFT calculations, it was then possible 
to determine the Hubbard parameters associated with the localized orbitals. For example, 
evaluating the energy of the state where there is one hole in the localized Cu (1^2 _y2 orbital 
(^Big [CuOq]~^'^) and using the state where there are no holes at all C'Aig [CuOq]~^^) 
as a reference yields the orbital energy for Cu dx2-y2. Evaluating the energy of the state 
where there are two holes in the Cu d^2_y2 orbital i^^ig [CuO^~^) then leads to the self- 
Coulomb term for this orbital. Similarly, evaluating the energy of a state where there is 
one hole in the plus combination of the localized Cu d^2_y2 and O p^j orbitals as compared 
to the state where there is one hole in the minus combination of these orbitals leads to the 
matrix element coupling the two orbitals. Clearly, all principle nearest neighbor Hubbard 
parameters describing the set of Cu, O, and apical O orbitals can be obtained in such a 
fashion using just this single CuOq cluster. 

The density that was optimized was for the undoped ground state, ^-Bi^ [CuOq]^^^ . 
This state has one hole in an orbital that is about 50% Cu dx2^y2 and 50% O Pcr- The state 
was chosen because it represents the closest approximation to the true density that can be 
obtained with this finite cluster. Note however that the point charge array reflects a doped 
state {x — 0.15) so the total charge on the system (cluster + point charges) is -0.3. While it 
was not possible to treat the cluster with a fractional charge to make the total charge on the 
system neutral, this discrepancy was effectively removed by the procedure described below 
in section II.B.l. 

The orbitals were localized using the Pipek-Mezey^'' localization procedure which max- 
imizes the sum of the squares of the atomic MuUiken populations over basis functions. This 
procedure was done in several steps. First, orbitals within a given irreducible representa- 
tion were localized and identified. This produced localized Cu d.j.2_y2 {Big) and dz2 {Aiy) 
orbitals, and symmetry combinations of the O (^ig, E.^, and Big), O p^^ {B2g, E.^, and 
A2g), and apical O Pz (^ig and A2u) orbitals. The procedure was then used on the symmetry 
combinations of the oxygen orbitals to obtain completely localized O Pa, O p^^, and apical 
O Pz orbitals. 

Hubbard parameters were derived by evaluating the DFT energies with fixed orbitals 
{i.e. non-SCF) as follows: 

1) The energy of the ^Aig [CuOq]~^^ state (having no holes) was evaluated. This was 

used as our reference state. 

2) The energies of the doublet [CmOs]^^^ states (each having a single hole in one of the 
localized orbitals) were evaluated. This yielded the orbital energies, Ef. 
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3) The energies of the singlet [CmOq] states (each having two holes in one of the localized 
orbitals) were evaluated. This yielded the self-Coulomb repulsion energies, Ja. 

4) The Hartree-Fock energies of the triplet and open-shell singlet [CuOq\^^ states (each 
having two holes in different localized orbitals) were evaluated. This yielded the exact 
exchange energy between orbitals, Kij. 

5) The energies of the triplet [CwOg]^^ states (each having two holes in different localized 
orbitals) were evaluated. This yielded the Coulomb repulsion energy between orbitals, 

Jy. 

6) The energies of the doublet [CuO^~^^ states (each having one hole in either the plus 
or minus combination of pairs of localized orbitals) were evaluated. This yielded the 
matrix elements coupling pairs of orbitals, Tjj. Note that the fully symmetric combi- 
nations of the oxygen orbitals were used in the evaluation of these terms rather than 
the completely localized orbitals. 

Parameters obtained from this procedure are listed in Table II. This parameter set will be 
referred to as the unsealed set. 

The Hubbard model that was developed included explicitly the Cu d^^-yi and dg2 
orbitals, the O orbitals (two per unit cell), and the apical O orbitals (two per unit 
cell). This led to a total of six bands. Solving the Hubbard model was done in two parts. 
The first took the input parameters (orbital energies and coupling terms) and found the 
single electron energies e„(A;) and wavefunctions (pn{k) where n is the band index. These k 
states were filled to the required doping level and the orbital occupations were evaluated. In 
the second step, these orbital occupations were used to reevaluate the orbital energies. The 
procedure was repeated until self-consistency was achieved. The Hubbard model is explained 
in more detail in the accompanying article.^ 

The Hubbard model was solved assuming no dispersion in the z direction. In calculating 
the density of states, however, z — axis dispersion was included perturbatively by assuming 
a dominant coupling through the apical O Pz orbitals. Further details follow. 

II. B. Orbital Energy Evaluation 
II.B.l. Coulomb Potential 

When orbital occupations are evaluated in the diagonalization step, the change in the 
Coulomb field had to be incorporated into the orbital energies, Ei. This is done by dividing 
up this field into that due to the CuOq cluster and that due to the point charge array. Since 
the reference orbital energies, E^, were originally defined for the case in which all CuOq 
orbitals were doubly occupied C'Aig [CuOq]^^^), the change in the orbital energies (£^j) due 
to the change in the occupation of the CuOq orbitals (iVj) were determined from, 

Ei^E^-{2- Ni)Jii - J2{2 - Nj){Jij - ^Kij), TV- > 1 (1) 
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Ei = Ef - Ju - ^(2 - Nj){Jij - -Kij), Ni<l (2) 

where the Cu dx^-y2 and (1^2 orbitals, the four O Pa orbitals and the two apical O Pz orbitals 
of the CuOq cluster were included in the summation. 

The Coulomb potential due to the changing long range field was evaluated by first 
subtracting off the potential due to the 1364 ion point charge array used in the DFT calcu- 
lations. This was evaluated as a classical point charge Coulomb interaction at each of the 
Cu, O, and apical O sites of the CuOq cluster. A similar Coulomb interaction was evaluated 
with a larger array having 22374 ions (17 unit cells wide in the a and b directions (64.396 A) 
and 5 unit cells tall in the c direction (66.030 A)). This was done to improve the long range 
Coulomb field over that which was used in the cluster calculations. This new Coulomb field 
was also broken up into components due to the Cu, O, apical O, and La/Sr sites. These 
fields were then appropriately scaled based on the orbital occupations from the Hubbard 
model and the effect was incorporated into the new orbital energies. 

II.B.2. Correlation 

An important point should be made here about the effect of correlation. As can be seen 
in equations (1-2), the self Coulomb interaction (Jjj) is treated separately from the other 
Coulomb interactions. This is a deviation from the mean field approximation of Hartree-Fock 
theory and conventional LDA band structure calculations. The mean field approximation 
would instead use the orbital energy correction equation: 

Ei^Ef-Y^(2-Nj)(Jij-^Kij) (3) 
j 

While equation (3) may be adequate for many materials, it breaks down in the limit of 
weakly interacting particles where the self Coulomb terms are much larger than the coupling 
matrix elements {Ja >> T^j). As can be seen from the data in Table II, this is the case 
for Lai g^SvQ i^CuO^. In this regime, equations (1-2) become valid. The difference between 
equations (1-2) and (3) can be seen when one orbital is at half occupancy {Ni — 1). In 
equation (2) the orbital energy is lowered by the full Ja term while in equation (3) the 
orbital energy is only lowered by (J^ — ^Ka) or equivalently ^Ju- Equations (1-2) assumes 
that correlation localizes all spins, while equation (3) inappropriately assumes the system has 
ionic character. With equation (3) there is a tendency to completely empty a band before 
electrons are removed from a second band in order to minimize the self-Coulomb energy. 
With the proper treatment of Ja in equations (1-2) there is instead a tendency to remove 
electrons from multiple bands in order to minimize the Coulomb repulsions between different 
orbitals. 

Variations on the type of correlation expressed in equations (1-2) have been introduced 
by several authors^^'^^ in studies of La2Cu04. Application of the self-interaction correction, 
or equivalently spin-polarization, has been successful in describing the antiferromagnetic 



5 



state of this material. As will be seen in this work, the introduction of correlation through 
equations (1-2) is critical to obtaining a band crossing in Lai ssSrQ ir^CuO^. 

II. B. 3. Orbital Relaxation 

Perhaps the main limitation in the method used here to obtain the Hubbard parameters 
is that the orbitals were not allowed to relax for different states. This will result in Coulomb 
interactions which are too high and orbital energies which arc too low. To determine the 
effect that orbital relaxation would have on our Hubbard parameters, we looked at isolated 
Cu, O, and apical O atoms in the full point charge array. We derived orbital energies and 
Coulomb interactions for these lone atoms using both fixed orbitals and fully optimized 
orbitals. The Coulomb energies were found to uniformly scale as 0.7 when relaxation was 
introduced. The orbital energies (defined here as E{(f?) — E{(f)) — Jjj) did not scale quite as 
uniformly, so different scales were used for Cu, O, and apical O. We found that the orbital 
energy for Cu scaled as 0.6, for O as 0.8, and for apical O as 0.7. We applied all these 
scales to the Hubbard parameters hsted in Table H to produce the set hsted in Table HI. 
It is this corrected set that was used in our final calculations. We should point out that 
we experimented with some other scaled sets (such as that obtained from uniformly scaling 
both the orbital energies and Coulomb terms by 0.7) and found no qualitative changes in 
the band structure. Even when the unsealed set was used, the same basic features were still 
observed. 

III. Results. 

In Figure 2 we show the 2-D dispersion of the top two bands as obtained with our 
Hubbard model using the scaled parameter set of Table HI. The band structure was computed 
for optimally doped La2-xSrxCu04^, where x — 0.15. This doping level corresponds to the 
removal of a total of 1.15 electrons per unit cell from the Cii/O/apical O bands (undoped 
La2CuO^ has 1 electron per unit cell removed from these bands). The orbital energies and 
occupations that were computed self-consistcntly from the model are given in Table IV. 

As can be seen in Figure 2, two bands appear to be important near the Fermi level. 
While no holes have been created in the lower band the Fermi level is just 0.035 eV above 
the top of this band. More importantly, a rigorous crossing, which is critical to the proposed 
theory, is seen just 0.153 eV below the Fermi level along the (0,0) — (tt/o, tt/o) symmetry 
line (note, there are actually four crossing points in the full 2-D Brillouin zone at {k,k), 
{k, —k), {—k, k), and {—k, —k)). The close proximity of this crossing to the Fermi level, as 
determined from this simple model, is a significant finding. Given errors in the Hubbard 
parameters and missing electronic effects, it is not difficult to conclude at this point that a 
band crossing, required of the IBP model, indeed appears to occur at the Fermi level of the 
optimally doped superconductor. This crossing can be characterized as arising between a 
broad "-Big" band dominated by Cu dy.2_y2 character and a narrower "Ai^" band dominated 
by Cu character. 

As shown in Figure 3, it is clearly the case that this band crossing is only observed 
when the self-Coulomb term is treated properly. Using the mean-field equation (3) we find 
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only a single isolated band having ''''Big" Cu d^^^yi character at the Fermi level. Other 
bands are well buried. We clarify this by noting that ligand field theory properly predicts 
that the Cu dx2-y2 orbital is the most unstable and d^^ is the next most unstable in the 
Jahn- Teller distorted octahedron of CuOq. The highest energy band orbitals of our Hubbard 
model reflect that this is true. However, as electrons are removed from the dj.2_y2 orbital the 
energy of this orbital stabilizes with respect to the dg2 orbital as a result of reduction in the 
self-Coulomb energy. This makes it more favorable to remove electrons from the d^2 orbital. 
Conversely, as electrons are removed from the (i^2 orbital the energy of this orbital stabilizes 
with respect to the d^^^y^ orbital. Yet because the self-Coulomb energy of the d^^ orbital is 
smaller than that of the dx2-y2 orbital due to s — (i^iz hybridization, it is actually possible to 
remove more electrons from the d.^2 orbital than the dj.2_y2 orbital. Such considerations are 
important since they lead directly to the observed band crossing. 

To take our calculations one step further, we considered what electronic effects might 
be missing that would have the greatest effect on the position of the crossing point relative 
to the Fermi level. Besides errors that might be present in the basic Hubbard parameters 
due to basis set hmitations, cluster size, and the DFT method itself, a number of missing 
key electronic effects can be identified. All are expected to lead to only minor perturbations 
of the band structure, but their cumulative effect could have an impact on the position of 
the band crossing. These effects include: 

1) Explicit inclusion of z-axis dispersion instead of the perturbative approach taken here. 

2) Inclusion of additional Hubbard parameters. Additional parameters should all be < 
0.04 eV. 

3) Explicit inclusion of other bands. Small mixings with other bands, in particular the 
Cu dxy/O Pn band, could have an effect on the position of the crossing. 

4) Use of more realistic charges in the cluster calculations. 

5) Inclusion of additional spin correlation. The current model only includes the self- 
interaction correction to account for the tendency of weakly interacting systems to 
localize spin. However, spin couplings between different orbitals, in particular the 
triplet coupling between Cu dx2-y2 and 0^22, has been ignored. 

6) Explicit inclusion of the La and Sr ions. The current model treats these atoms as 
point charges having the formal charges -1-3.0 and -1-2.0 , respectively. In reality, these 
ions are hkely less highly charged. They also have spatial extent which leads to Pauli 
repulsions. 

7) Inclusion of the orthorhombic distortion. The crystal in the superconducting phase 
is distorted from its high symmetry tetragonal (-D4/1) structure to a lower symmetry 
orthorhombic (C2/1) structure.^ This is manifested by a tilting of the CuOq units or 
equivalently a buckling of the Cu02 planes. An important effect from this distortion is 

that two of the four crossing points (at {k, —k) and {—k, k)) become strictly avoided due 
to the reduced symmetry. However, it should be stressed that the other two crossing 
points (at {k,k) and {—k, —k)) are maintained, which is critical to the proposed IBP 
theory. 
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While it is not clear that all of these effects will be favorable in terms of moving the band 
crossing toward the Fermi level, their combined effect could easily lead to such a change. 
We stress that a perturbation of only 0.153 eV (a small quantity on the chemical scale) is 
necessary to observe a Fermi level crossing. 

We have, in fact, applied a specific perturbation to our model in order to incorporate 
one of the above effects. Investigations of the effect of La and Sr using a variety of clusters 
at the BLYP level showed significant mixing of the O orbitals with the orbitals of these 
metals. The implication of these results is that there is diffusion of the O p-j^ electrons onto 
the La/Sr sites which has been ignored by the CuOq cluster calculations. The effect of 
this charge diffusion is to lower the energy of both the O p^^ and O p^^ orbitals and raise 
the energy of the apical O Pz orbitals relative to the Cu orbitals. To account for this effect 
we included in our Hubbard model calculations an adjustable parameter which defined the 
extent of charge transfer from O pj^ to La/Sr. While we do not explicitly include these tt 
bands in the final model, Coulomb terms for the O orbitals were evaluated and included 
in Tables II and III (La and Sr were still treated as classical point charges). This allowed 
us to include the effect of charge transfer by altering the Coulomb field in a fashion similar 
to that explained in section II.B.l. Orbital energies were reevaluated to reflect the change 
in the Coulomb held due to this charge transfer. 

The primary effect of this charge transfer is to stabilize the tt bands relative to the 
Fermi level by significantly lowering the orbital energy oi O p-j^. A secondary effect, however, 
is to raise the crossing point of the two bands of interest closer to the Fermi level for the op- 
timally doped system. This was accomplished when the charge transfer term was empirically 
adjusted to the value of 0.50 electrons transfered. That is, the charge on the La/Sr sites 
was +2.425 compared to the formal charge of +2.925. This should be considered within the 
range of reasonable charges for these ions (note the charges on the other atoms are +1.563 for 
Cm, -1.429 for O, and -1.777 for apical O). We should note though, that the extent of charge 
transfer calculated here may be an overestimate in light of the fact that the other electronic 
effects hsted above have not yet been incorporated. But, as detailed in the accompanying 
article by Tahir-Kheh,^ the resulting band structure proves to have the necessary features 
to explain a number of key experiments. 

The 2-D dispersion of the top two bands from these calculations is shown in Figure 
4 and the optimized orbital energies and occupations are given in Table V. In addition the 
Fermi surfaces are shown in Figure 5 and the density of states of the top two bands is shown 
in Figure 6. 

IV. Discussion. 

IV. A Band Structure of Lai.85Sro.i5Cu04 

The band crossing we observe arises from the following considerations. In this discus- 
sion, we will refer to the two bands that cross as the "-Big" and "^ic," bands. The crossing 
produces two new bands that touch, referred to as t/(pper) and L(ower). The dispersion of 
the "-Big" band, dominated by Cu dx'^_y2 character, is rather broad, on the order of 2 eV, 
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producing a low density of states. The dispersion of the "^ig" band, dominated by Cu d^^ 
character, is in contrast rather narrow, on the order of 0.3 eV, producing a high density of 
states. At the (vr/a, vr/a) point, the higher energy band is ''Big" in nature, corresponding to 
the completely antibonding combination of Cu dx2-.y2 and O orbitals. The lower energy 
band is "^i^" in nature, corresponding to the antibonding combination of the Cu dg2, O p^, 
and apical O orbitals. At the (0,0) point, however, the higher energy band is "^i^" in 
nature and the lower energy band is "-Big. Due to this change in the relative energetics 
between the top of both bands and the bottom of both bands, the "-Bi^" and "^ig" bands 
must cross, producing a U band and L band which each have both "Big" and "^ig" char- 
acter. This crossing of bands is strictly avoided everywhere except for one point in the 2-D 
Brillouin zone. Along the (0, 0) — (tt/o, tt/o) symmetry line, the crossing is rigorously allowed 
since the ''Big" and "Aig" orbitals cannot mix along this direction. Due to this rigorous 
crossing, the U band and L band must touch. While the position of this crossing point (or 
touching point) is subject to variation, its existence is quite robust over a wide range of 
model parameters. 

In the closeup of the density of states shown in Figure 6, it can be seen how the crossing 
of the "Big" and "Aig" bands affects the nature of the U and L bands. The U band starts 
at +0.47 eV and represents the (tt/o, tt/o) point of the "Big" band. The density of states of 
this band remains consistently low until about +0.10 eV where there is a sharp peak. This 
peak represents the change in character of the band from "Big" to "Aig" in the vicinity of 
the (7r/a,0) point. The density of states of the U band remains relatively high at energies 
below this point, being dominated by "Aig" character. The band terminates at -0.25 eV, 
which represents the (0,0) point of the "Aig band. The onset of the L band is characterized 
by a sharp peak in the density of states at +0.03 eV. This peak occurs in the vicinity of the 
{n/a,7i/a) point of the "Aig" band. The peak comes down to a low density of states near 
the Fermi level and the density of states remains low, being dominated by "Big" character 
until about -0.20 eV when character from other orbitals (in particular the "A2u" antibonding 
combination of the apical O p^ orbitals) starts to mix in. The change in character of this 
band at lower energies should have no effect on the issue of the band crossing or the IBP 
model. 

The character of these bands can be seen more clearly in Figure 7. In Figure 7a, the U 
band is shown to be dominated by Cu d^^i-y^ character above 0.10 eV while the L band is 
dominated by d,j.2_y2 below the Fermi level. In Figure 7b, the U is shown to be dominated 
by Cu dz2 character below 0.10 eV while the L band is dominated by dg2 only at the Fermi 
level. 

While most band structure calculations have only shown a single band at the Fermi 
level,^^ we argue that this is due to the lack of correlation. In fact, the work of Shiraishi, et 
al.,^^ which included the effect of spin-polarization, found that doping of La2Cu04 resulted 
in the formation of two types of holes from two distinct bands. While no band crossing was 
noted, the two types of holes were characterized as being Cu drj.2_y2 and rf^z. Eto, et al}"^ 
also suggested the importance of the Cu d^^ orbital based on cluster model calculations. We 
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find that this work lends support to the findings reported here. 

IV. B. Band Structure of Related Materials 

Certain qualitative features in the density of states for La2-xSrxCu04 (LASCO) 
should be characteristic of many superconducting copper oxide materials, such as 

YBa2Cu30Q+s (YBCO), Bi2Sr2Can-iCun02n+6+y, Tl2Ba2Can-iCun02n+'i+5, and oth- 
ers. Each should be characterized by a peak in the U band just above the top of the L band. 
This is due to band repulsions away from the (0,0) — (n/a.-ri/a) diagonal which introduce 
"y4ip" character into the U band. When these band repulsions are large enough (such as 
they are at (vr/a, 0)), the peak in the density of states of the U band should occur above the 
top of the L band. The L band, on the other hand, should have a sharp peak at its onset at 
(7r/a,7r/a) which vanishes as the band becomes "Big" in character. 

Although the electronic structure of YBCO should share common features with 
LASCO, it differs in one important regard: the local symmetry of YBCO in the super- 
conducting state is -D2/1 as compared to C2h for LASCO.^ This symmetry will not allow a 
rigorous crossing of the two bands along the diagonal (or any other point) in the 2-D Bril- 
louin zone. However, it should be recognized that the dual Cu02 planes of YBCO lead to 
four bands in contrast to the two bands of LASCO. Dispersion of these four bands in the z 
direction leads to two bands having "A^" symmetry at A;^ = and kg = n/c and two bands 
having "Bi^" symmetry. The two "^g" bands are principally composed of the bonding com- 
binations of the Cu d^2_y2 orbitals from the two planes and the bonding combinations of 
the Cu d^2 orbitals from the two planes. The bands are the antibonding analogues. 

While the two "Ag" bands are not precluded from mixing at any symmetry point, and the 
two "Biu" bands are also not precluded from mixing, the "Ag" and "-Bi„" bands cannot mix 
when = or kz = tt/c. This suggests that crossings could occur between the "Ag" bands 
and the "-Biu" bands. Based on this idea, we propose the following scenario. (Note, we do 
not address the effects of the chain CuO bands. These bands will perturb the planar CUO2 
bands, but we suspect the essential topology described below holds.) 

There are two sets of two bands {U and L) having qualities similar to those shown 
here for LASCO. One set has "A^" symmetry while the other has "Si^" symmetry. The 
primary difference between these bands and those of LASCO is that the crossing along the 
(0,0) — (7r/a,7r/a) direction is avoided. Since the Cu d,j.2_y2 dispersion in the z direction is 
expected to be small, the d^2_y2 components of these bands (analogous to the bands 
of LASCO) should be nearly degenerate. In contrast, the Cu dz2 components of these bands 
(analogous to the "Ai^" bands of LASCO) should be separated in energy. 

The onset of superconductivity in YBCO should occur when the Fermi level coincides 
with a band crossing. This hkely occurs a,t k^ — n/c between the "-Bi^" L band (which is 
completely antibonding in the z direction) and the "Ag" U band. Interestingly, it is not 
required to occur along the (0,0) — (7r/a,7r/a) direction. In fact, there is likely a double 
crossing within a quadrant of the Brillouin zone, as depicted in Figure 8. In this schematic, 
we view the Fermi surfaces as squared circles centered around (0, 0) for "-Bi^" L and "A^" 
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L and centered around {n/a^Ti/a) for "-Bi^" U. Superconductivity begins when the "A^" L 
and "-Biu" U bands first touch. Since the radius of the "-Bi^" U Fermi surface is increasing 
faster than the radius of the "Ap" L Fermi surface is decreasing (due to the difference in 
their densities of states), two crossing points can be sustained over a wide doping range. At 
sufficiently higher doping levels, a second set of band crossings of a similar nature should 
occur at — 0. This proposal easily explains the extended doping range observed for YBCO 
and the appearance of two T^'s in different doping regimes for this material. Furthermore, 
it suggests a reason for the increasing for LASCO vs. underdoped YBCO vs. optimally 
doped YBCO. We correlate an increase in the number of crossing points, or, more correctly, 
the number of crossing points which are thermally accessible, to an increase in T^. Indeed, 
it is observed that superconductivity begins with YBa2CuzOQ_Q showing a of 60 K. This 
Tc is sustained upon further doping until a rapid increase to Tc = 90 K is observed near 
YBa,CuA,. Farther doping to YBa,Cu,0,l maintains at thte higher temperature.'^ 

Similar analysis can be applied to the bismuth^^ and thallium^^ systems. Wc note that 
in these systems an increase in is correlated with an increase in the number of Cu02 
planes per unit cell. Following the above arguments this makes sense in that it leads to more 
bands which produce more band crossings as shown schematically in Figure 9. For a three 
plane system we are now dealing with six bands. Four of these bands will be symmetric with 
respect to reflection through the middle plane (bands gl, g2, g3, and g4) and two of these 
bands will be antisymmetric (bands ul and u2). We argue that band repulsions near the 
(vr/a, 0) point, as seen in the LASCO band structure, would be stronger for the g bands, 
since dg2 character should appear in these bands at a higher energy than in the u bands. 
This could produce a crossing between the ul and g2 bands at A;^ = and kz — tt/c as 
depicted in the figure. Since the two surfaces might be expected to be highly coincident at 
this point, the number of thermally accessible crossing points should be high. Clearly, it can 
be seen that the addition of more Cu02 planes increases the probability of favorable crossing 
situations as illustrated here, and this should lead to potentially higher T^'s. 

Finally we wish to note that the electron doped system, Nd2-xCexCu04,^^ may be 
substantially different in regard to the nature of the two bands that cross. While we an- 
ticipate one of the bands will be "-Bi^" Cu d^2_y2 in character, the other band is likely 
not "^ig" Cu d^2 since this band only appears as electrons are removed from the system. 
We suggest instead that the second band is a Nd/Ce band. It follows from this suggestion 
that high temperature superconductivity is not dependent on the specific crossing of the 
Cu dx2-y2 and dg2 bands. We stress, in fact, the only requirement of the IBP model that 
we now see is that some sort of crossing of bands occurs at the Fermi level. This leads us 
to be optimistic that with the careful exploitation of symmetry, entirely new classes of high 
temperature superconductors will be developed in our future. 

V. Conclusion. 

We have presented the results of a new Hubbard model calculation on the optimally 
doped superconducting material Lai g^SrQ ir^CuO^. We conclude from these calculations 
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that there is a crossing of two bands which occurs at the Fermi leveL One of these bands is 
''Big" in character, dominated by Cu dx2_y2 and the other is "^ig" in character, dominated 
by Cu d^2. The crossing rigorously occurs along the (0, 0) — (vr/a, vr/a) symmetry line of the 
2-D Brillouin zone. As detailed in the accompanying article by Tahir-Kheli,^ an inter-band 
pairing (IBP) of electrons between these two bands leads to superconductivity. It can only 
occur at the critical doping level where the Fermi energy coincides with the band crossing, 
i.e. X = 0.15. The resulting density of states from this work is used in the accompanying 
work to explain a number of key experiments on this material. 

Extension of the band model obtained for LASCO in these calculations to Y BCO 
provides an easy explanation for the observation of a wide doping range for superconductivity 
to occur in this material as well as an explanation for the observation of two T^'s in different 
doping regimes. We anticipate at this point that other materials will be similarly understood 
largely through geometric considerations. 

Acknowledgment: The authors wish to thank Dr. Jean-Marc Langlois for many useful 
discussions. 
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Table I. Crystal structure of Lai^s^SrQ^i^CuO^ (in A). 



a 


3.788 


b 


3.788 


c 


13.206 


Cu 


(0.000,0.000,0.000) 


0(1) 


(0.500,0.000,0.000) 


0(2) 


(0.000,0.000,0.182) 


La/Sr 


(0.000,0.000,0.361) 



Table II. Unsealed Hubbard parameters (in cV). E is an orbital energy, T an orbital coupling 
matrix element, J a Coulomb repulsion term, and K an exchange energy term. 



E(x2-y2) 


-3.109 




25.563 


K{x^-y^/z'') 


1.219 




-3.338 


J(a;2 - 1/2/0 p,) 


7.999 


Kix^ - y^O p,) 


0.081 


E{0 p,) 


-10.417 


J(x2 - y2/02 p,) 


6.471 


K(x2 - y2/02 p,) 


0.028 


E{02 p,) 


-12.442 


J{x^-y^/Op,) 


7.053 


K{x^ - y^O p,) 


0.006 


J{x^ - y^/x^ - y^) 


29.281 


Jiz^Opa) 


6.471 


K{z^lO Pa) 


0.028 




25.990 


J(^2/02 p,) 


6.855 


K{z^l02 p,) 


0.262 


J{0 pJO p,) 


17.375 


J{Z^I0 Pn) 


6.813 


Kiz^/0 p^) 


0.008 


J(02 p,/02 p,) 


11.365 


J(0 pJO p',) 


5.345 


K{0 pJO p',) 


0.035 


T{x^-y^/Op,) 


1.347 


AO pjo p'^) 


3.890 


K{0 pJO p'^) 


0.008 


T{z^/0 p,) 


0.514 


J{0 pJ02 p,) 


5.144 


K{0 pJ02 p,) 


0.211 


r(zV02 p,) 


1.076 


J{0 Pa/Op,) 


14.299 


K{0 pjOp,) 


0.677 


no pJO p',) 


0.368 


J{0 pJO p'^) 


5.548 


K{0 pJO p'^) 


0.092 


no pJO p'^ 


-0.041 


J{0 pJO p'i) 


3.953 


K{0 pJO p'i) 


0.008 


no pjo2 p,) 


0.078 


J(02 pJ02 p',) 


3.502 


K{02 pJ02 pi) 


0.098 


T(02 p,/02 p',) 


0.493 


J(02 p,/0 p,) 


4.868 


K{02 p,/0 p,) 


0.087 
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Table III. Scaled Hubbard parameters (in eV). E is an orbital energy, T an orbital coupling 
matrix element, J a Coulomb repulsion term, and K an exchange energy term. 



r^yx y ) 


1 .uuo 


J [X y / z ) 


1 7 RQA 


j\[x y / z ) 


1 91 Q 


J^{Z ) 


n ^Qf, 


•J\X U 1^ Pa) 


O.OcJcJ 


KiT^ 7/2 in n \ 

j\\x y 1^ Pa) 


U. UO -L 


Pa) 


-lU.U < 1 


'jy^ ~y 1^^ Pz) 




~y Pz) 


U.UZo 


tj\Ul Pz) 


-o. / uy 


7/'^2 .,2 //-) ^ \ 

j[x -y /u PttJ 


4.yo/ 


K[x -y /U PttJ 


U.UUD 


J{x -y /x -y ) 


on A HT 

zU.4y / 


J[Z /Opa) 


o 1 o 
o.olo 




O.Ozo 




18.193 


J{zV02 Pz) 


4.799 


Kizy02 p,) 


0.262 


J{0 pJO p,) 


12.163 


J{z''/0 p,) 


4.769 


Kiz^O p^) 


0.008 


J{02 pJ02 Pz) 


7.956 


J{0 pJO p',) 


3.742 


K{0 Pa/0 p',) 


0.035 


T{x^ - y^O p,) 


1.347 


J{0 pJO p'^) 


3.601 


K{0 Pa/0 p'/,) 


0.008 


T{zyo p,) 


0.514 


J{0 p,/02 Pz) 


3.601 


K{0 Pa/02 Pz) 


0.211 


T{z^l02 p,) 


1.076 


J{0 pjO p^) 


10.009 


K{0 Pa/ Op,) 


0.677 


no pjo p',) 


0.368 


AO Pa/0 p'^) 


3.884 


K{Opa/0 p'^) 


0.092 


no pjo p'!,) 


-0.041 


J{0 pJO p'i) 


2.767 


K{o Pa/0 p';) 


0.008 


no Poioi Pz) 


0.078 


J{02 pJ02 p'z) 


2.451 


K{02 Pz/02 p'z) 


0.098 


T{02 Pz/02 p'z) 


0.493 


J{02 Pz/0 p,) 


3.408 


K{02 Pz/0 p^) 


0.087 



Table IV. Computed orbital energies and orbital occupations for optimally doped 
Lai,s5SrQ ic,Cu04 taken from our Hubbard model using the scaled parameter set of Ta- 
ble HI and no O p^^ to La/ Sr charge transfer. All energies are relative to the Fermi level (in 
eV). 



Orbital 


Energy 


Occupat: 


Cu x^ — y^ 


-2.570 


1.572 


Cu z^ 


-1.663 


1.785 


Pa 


-4.254 


1.806 


02 Pz 


-1.611 


1.941 
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Table V. Computed orbital energies and orbital occupations for optimally doped 
Lcii.85SrQ i^Cu04 taken from our Hubbard model using the scaled parameter set of Ta- 
ble III and 0.50 electron O to La/Sr charge transfer. All energies are relative to the 
Fermi level (in eV). 



Orbital Energy Occupation 

Cux'^-y'^ -2.403 1.770 

Cu -2.092 1.666 

Opa -6.122 1.929 

02 p^ -0.852 1.777 
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FIG. 1. Tetragonal unit cell of Lai.85Sro.i5Cu04 show- 
ing the CuOe cluster used in the DFT calculations in bold. 




(0,0) 



(7r/a,7r/a) (7r/a,0) 



(0,0) 



FIG. 2. Dispersion of the top two bands obtained 
from the Hubbard model using the scaled parameter set of 
Table III and no O Ptt to La/Sr charge transfer, as shown 
along the symmetry line (0,0) — (7r/a,7r/a) — (7r/a,0) — 
(0,0). The Fermi level is shown to lie just above the top 
of the lower band and only 0.153 eV above the crossing 
point along the (0,0) — (77/0,77/0) symmetry line. 



1.5 




(0,0) (7r/a,7r/a) (7r/a,0) (0,0) 



FIG. 3. Dispersion of the top band when the mean- 
field equation (3) is used in the Hubbard model (see de- 
scription of Figure 2). Other bands are several eV lower 
in energy and are not shown. 




(0,0) (Tr/a,0) (7r/a,7r/a) (0,0) 



FIG. 4. Dispersion of the top two bands obtained 
from the Hubbard model using the scaled parameter set of 
Table III and 0.5 electron O Ptt to La/Sr charge transfer, 
as shown along the symmetry lines (0,0) — (7r/a,7r/a) — 
(tt/o, 0) — (0, 0). The Fermi level is shown to coincide with 
the crossing point along the (0,0) — (7r/o,7r/a) symmetry 
line. 



7r/2 







7r/2 



FIG. 5. Fermi surfaces showing the touching of the 
Upper and Lower bands along the (0. 0) — (tt/o, nja) sym- 
metry line. The squareness of this surface is expected to 
change with the addition of longer range terms. 
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FIG. 6. Density of states of the top two bands shown 
from E = -0.25 eV to B = +0.50 eV. 
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FIG. 7a. Extent of Cu d^i_y2 orbital character in the 
Upper (solid line) and Lower (dashed line) bands. 
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FIG. 7b. Extent of Cu d^-i orbital character in the 
Upper (solid line) and Lower (dashed line) bands. 
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FIG. 8. Schematic depiction of Fermi surfaces for YBa2Cu306+5. The onset 
of superconductivity is expected when the "Bi„" L band and the "Ag" L band 
first cross at hz = 'Kjc (top). Superconductivity continues upon further doping 
(bottom). A change in Tc may be expected when a similar set of crossing appears 
for kz = 0. The YBCO unit cell is shown (left). 




FIG. 9. Schematic depiction of Fermi surfaces for a typical 3-plane copper 
oxide superconductor. Band repulsions may lead to highly coincident g2 and u\ 
surfaces in the vicinity of {T:ja,Q) when kz = d ov kz = i:lc. 



